Nuclear-spin dependent parity violation in diatomic molecular ions 
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Nuclear-spin-dependent (NSD) parity violating (PV) effects can be strongly enhanced in diatomic 
molecules containing heavy atoms. Future measurements are anticipated to provide nuclear anapole 
moments and strength constants for PV nuclear forces. In light molecules, the NSD electroweak 
electron-nucleus interaction may also be detected. Here we calculate NSD PV effects for molecular 
ions. Our calculations are motivated by rapid developments in trapping techniques for such systems 
at low temperatures. 
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It was previously shown that nuclear spin-dependent 
(NSD) parity violation (PV) effects are enhanced by a 
factor of 10 5 in diatomic molecules with 2 E 1/ / 2 and 2 Hi/2 
electronic states due to the mixing of close rotational 
states of opposite parity (fi-doublet for 2 n 1 / 2 ) [IH3]- De- 
Mille and co-workers suggested measuring NSD PV ef- 
fects by using neutral diatomic molecules in a Stark in- 
terference experiment to determine the mixing between 
opposite-parity rotational/hyperfme levels [1 . Another 
proposal was published in Ref. [S], and correspond- 
ing experiments have already started. Recently, it was 
demonstrated that positive molecular ions may be easily 
trapped and studied at low temperatures S6 ], which mo- 
tivated us to perform calculations on NSD PV effects in 
such systems. 

The term in the Hamiltonian operator arising from the 
NSD parity violating electron-nucleus interaction is 

H A = k NS d^^y~p(t). (1) 

Here, and throughout the text, we use atomic units. In 
Eq. 0, knsd is the dimcnsionless strength constant, 
Gp = 2.22249 x 10~ 14 a.u. is the Fermi constant, a is a 
vector comprised of the conventional Dirac matrices, I is 
the nuclear spin, r is the displacement of the valence elec- 
tron from the nucleus, and p(r) is the (normalized) nu- 
clear density. There are three sources for this interaction: 
the first contribution arises from the electroweak neu- 
tral coupling between electron vector and nucleon axial- 
vector currents (V e Ajv) [?]■ The second contribution 
comes from the nuclear-spin-independent weak interac- 
tion combined with the hyperfine interaction [5] . Finally, 
the nuclear anapole moment contribution, which scales 
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with the number of nucleons A as ka ~ A 2 / 3 , becomes 
the dominant contribution in spin-dependent atomic PV 
effects for a sufficiently large nuclear charge Z [9j [10] . 

The anapole moment was first predicted by Zeldovich 
[TT] in 1958 as a new parity (P) violating and time (T) 
reversal conserving moment of an elementary particle. It 
appears in the second-order multipole expansion of the 
magnetic vector-potential simultaneously with the P— 
and T- violating magnetic quadrupole moment [12] . The 
nuclear anapole moment was experimentally discovered 
in the 133 Cs atom in 1997 [13] following a proposal by 
Flambaum and Khriplovich [H] , who showed that the nu- 
clear anapole provides the dominant contribution to the 
nuclear-spin-dependent parity violating effect in atoms 
and molecules. 

The nuclear anapole requires nuclear spin I ^ and 
in a simple valence model has the following value |10) . 

k a = 1.15 x 1(T 3 (7^7) A 2 ' 3 ^. (2) 

Here, K = + 1/2), I is the orbital angular 

momentum of the external unpaired nucleon i — n,p; 
Hp = +2.8, jtx„ = —1.9. Theoretical estimates give 
the strength constant for nucleon-nucleus weak potential 
g p m +4.5 for a proton and \g n \ ~ 1 for a neutron [14] . 
The aim of anapole measurements is to provide accurate 
values for these constants, thus obtaining important in- 
formation about hadronic weak coupling. 

A number of theoretical investigations of the nuclear 
spin-dependent parity violation in diatomic molecules 
have been performed in recent years, using both semiem- 
pirical [U irSHTT] and ab initio methods ITM2"2"] . In a 
recent paper |23] we presented Dirac Hartree-Fock and 
relativistic density-functional calculations of the elec- 
tronic Wa factor of the diatomic group-2 and -12 fluo- 
rides and a number of other diatomic compounds. In this 
work we investigate the nuclear spin-dependent parity vi- 
olation effects in a different type of system, i.e. positively 
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charged dimers, as these systems have an experimental 
advantage of being easily trapped [6 . Diatomic ions have 
also been proposed for the search for the electron electric- 
dipole moment (eEDM) j2"4"H2l)] , and the preliminary ex- 
periments are currently being conducted [27] ■ Here we 
use the combination of methods presented in Ref. [23] to 
calculate the Wa factors of positively ionized group-13 
and group-15 fluorides ( 2 S 1/ / 2 and 2 n 1 / 2 ground states, 
respectively), and a number of other positive diatomic 
ions having a 2 £i/2 ground state. 

For 2 £;l/2 and 2 n 1 / 2 electronic states, the interaction 
can be replaced by the effective operator, which ap- 
pears in the spin- rotational Hamiltonian [Hill], 
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(n x S') • I 
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where S' is the effective spin and n is the unit vector 
directed along the molecular axis from the heavier to the 
lighter nucleus. The electronic factor Wa is found from 
evaluating the matrix elements of the ctp(r) operator in 
the molecular spinor basis [55]. The 2 £i/2 and the 2 n 1 / 2 
open-shell electronic states are two-fold degenerate, cor- 
responding to the two possible projections of electronic 
angular momentum along n, i.e. \ — |±i). When oper- 
ating within this degenerate space, the operator ^ap(r) 

is equivalent to W/i(nX S') (Eq. Time-reversal sym- 
metry ensures that only the matrix elements that are off- 
diagonal in f2 are non-vanishing. This symmetry rule is 
encapsulated within the effective operator H A S by the 
angular factor (n x S'). Here the effective spin S' gener- 
ates rotations in the degenerate subspace analogously to 
usual spin operator S in a spin-f/2 system. 

The calculations were carried out within the open- 
shell single determinant average-of-configuration Dirac- 
Hartree-Fock approach (DHF) [55] and within the rela- 
tivistic density functional theory (DFT) [30 , employing 
quaternion symmetry |31[ 132] . A finite nucleus, mod- 
eled by the Gaussian charge distribution was used [55] , 
All the calculations were performed using the developer's 
version of the DIRAC10 program package [34] , 

For the lighter elements (boron to phosphorus), uncon- 
tracted aug-cc-pVTZ basis sets were used [551 ES] • For 
the rest of the atoms, we employed Faegri's dual fam- 
ily basis sets [37] . As a good description of the electronic 
wave function in the nuclear region is essential for obtain- 
ing reliable results for parity violating properties [38] , we 
augmented the basis sets with high exponent s and p 
functions, which brings about an increase of around 10% 
in the calculated values of Wa- The basis sets were in- 
creased, both in the core and in the valence regions, to 
convergence with respect to the calculated Wa constants. 
The final basis sets can be found in Tabic fl] 

Where available, experimentally determined bond dis- 
tances R e were used. For molecules where R e is not 
known experimentally, we optimized the bond distance 
using relativistic coupled cluster theory with single, dou- 
ble, and perturbative triple excitations, CCSD(T) [55] , 



TABLE I. Basis sets employed in the calculation of the Wa 
constants. All elements with Z > 15 are described by the 
Faegri basis sets [37] augmented by high exponent, diffuse, 
and high angular momentum functions. 



Atom 



Basis Set 



B 

N 

O 

F 

Al 

P 

Ga 

As 

Y 

Zr 

In 

Sb 

Hf 

Tl 

Bi 

Ac 



5 
7 

8 
9 

13 
15 
31 
33 
39 
40 
49 
51 
72 
81 
83 



aug-cc-PVT5 
aug-cc-PVTZ 
aug-cc-PVTZ" 
aug-cc-PVTZ 
aug-cc-PVTzp 
aug-cc-PVT^ 
22sl9pWd8f2g 
21s20pllrf8/2g 
21s20pl2rf9/2g 
21s20pl2rf9/2g 
22s20pl2rf9/2g 
22s21pl3rf9/2g 
25s22pl6rfl0/2ff 
25s23pl5rfl0/2ff 
25s24pl6rfll/2t/ 
26s24pl6<m/2p 



a augmented by 3 high exponent p functions 
b augmented by 4 high exponent s and 3 high exponent p 
functions 

c augmented by 4 high exponent p function. 
d augmented by 1 high exponent s and 3 high exponent p 
functions 



To reduce the computational effort, wc employed an in- 
finite order two-component relativistic Hamiltonian ob- 
tained after the Barysz-Sadlej-Snijders (BSS) transfor- 
mation of the Dirac Hamiltonian in a finite basis set 
[4U1 |4"T] . Our calculated R e are typically within 0.01 A 
of the experimental values, where available. The experi- 
mental/calculated equilibrium distances can be found in 
Table M 

In the DFT calculations we used the Coulomb- 
attenuated B3LYP functional (CAMB3LYP*), the pa- 
rameters of which were adjusted by Thierfelder et al. |45] 
to reproduce the PV energy shifts obtained using coupled 
cluster calculations (the newly adjusted parameters are 
a = 0.20, p = 0.12, and // = 0.90). 

In our previous work [23] we have examined and 
compared various schemes for adding electron correla- 
tion to the Dirac-Hartree-Fock Wa values, and core- 
polarization contributions to the DFT results. Here, we 
correct the calculated DHF and DFT Wa for core polar- 
ization using a scaling parameter, Kcp- This parameter 
is obtained from atomic calculations as described in the 
following. The main contribution to the matrix elements 
of the NSD interaction for the valence molecular elec- 
trons comes from short distances around the heavy nu- 
cleus, where the total molecular potential is spherically 
symmetric to very high precision, and the core of the 
heavy atom is practically unaffected by the presence of 
the second atom, justifying our use of the atomic model. 
The molecular orbitals of the valence electron can thus 
be expanded in this region, using spherical harmonics 
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TABLE II. Internuclear distances R e (taken from 
CCSD(T) calculations, unless referenced otherwise, A), 
core-polarization scaling parameters Kcp, the P-odd 
interaction constants Wa (Hz) obtained using DHF 
and DFT, and the final recommended values, taken as 
WA(Final)= {Wa(T>¥T)K CP + W A (DEF)K C p)/2. Rela- 
tivistic factors Rw (see Eq. 1 10}) are also shown. 



Z R w R e (A) K C p 



W a (Hz) 



BF+ 

A1F+ 

GaF+ 

InF+ 

T1F+ 



Group 13 ( 
5 1.01 1.314 1.1 
13 1.07 1.59C0 1.2 
31 1.41 1.68c 5 1.1 
49 2.2 1.91 1.1 
81 7.4 2.00 1.1 



DHF 


DFT 


Final 


^1/2) 






1.74 


1.71 


1.90 


9.62 


10.39 


12.0 


94.4 


93.5 


103.4 


370.3 


358.3 


400.7 


3833 


3622 


4100 



NF+ 


7 


1.02 1.18(f 




-0.014 


-0.015 


-0.016 


PF+ 


15 


1.10 1.524^ 


1.2 


-0.16 


-0.18 


-0.20 


AsF+ 


33 


1.47 1.660 


1.1 


-6.48 


-7.29 


-7.67 


SbF+ 


51 


2.3 1.832 


1.1 


-60.0 


-66.1 


-71.6 


BiF+ 


83 


8.1 2.281 


1.1 


-2204 


-2123 


-2380 






Other systems ( 2 Y} 1 / 2 ] 






YF+ 


39 


1.69 1.885 


1.2 


45.6 


39.4 


51.0 


ZrO+ 


40 


1.73 1.695 


1.2 


43.2 


31.1 


44.6 


HfO+ 


72 


5.0 1.708 


1.2 


662.5 


609.1 


762.9 


AcF+ 


89 


11 2.106 


1.2 


1654 


1614 


1961 



a Ref. 
b Ref. 
c Ref. 



42 
43 
44 



centered at the heavy nucleus, 

\ip v ) = a\s 1/2 ) + b\p 1/2 ) + c\p 3/2 ) 



d\d 3/2 



(4) 



Only S1/2 and P1/2 terms of this expansion give signif- 
icant contribution to the matrix elements of the weak 
interaction. These functions can be considered as states 
of an atomic valence electron and are calculated using 
standard atomic techniques in two different approxima- 
tions: one that includes electron correlation and another 
that does not. 

The single electron DHF Hamiltonian is given by 



H Q = ca • p + (f3 - iy -- + V e (r), 

r 



(5) 



where a and /3 are the Dirac matrices and V e {r) is the 
self-consistent DHF potential due to atomic electrons. 

The self-consistent DHF procedure is first performed 
for the closed shell ion, from which the valence electron is 
removed. Then the core potential V^^'" is frozen and 
the valence Sx/ 2 and pi/ 2 states are calculated by solving 
the DHF equation for the valence electron, 



(H - e v )ip v = 0, 
where Ho is given by ib\. 



(6) 



The core polarization can be understood as the change 
of the self-consistent DHF potential due to the effect of 
the extra term (the weak interaction operator Ha) in 
the Hamiltonian. The inclusion of the core polarization 
in a self-consistent way is equivalent to the random-phase 
approximation (RPA, see, e.g. [46]). The change in the 
DHF potential is found by solving the RPA-type equa- 
tions self-consistently for all states in the atomic core, 



(H -e c )5i> c = -(H A + 5V A )i; c 



(7) 



Here, Ho is the DHF Hamiltonian |5|), index c enumer- 
ates the states in the core, 5ip c is the correction to the 
core state c due to weak interaction Ha, and SVa is the 
correction to the self-consistent core potential due to the 
change of all core functions. Once SVa is found, the core 
polarization can be included into a matrix element for 
valence states v and w via the redefinition of the weak 
interaction Hamiltonian, 



(v\H A \w) ->■ (v\H A + SV A \w) 



(8) 



We then obtain the scaling parameter for core- 
polarization effects, Kqp, from 



A', 



CP 



{^ 2 \H A \^ 1/2 ) 



(9) 



It should be noted that for the positively charged group 
15 fluorides we have only calculated the correlations be- 
tween the valence electrons and the core; the correlations 
between the valence ns and np electrons are not included. 

We investigated two types of positively ionized di- 
atomic molecules: those with a 2 Si/2 ground state, in- 
cluding group 13 fluorides and a number of other sys- 
tems, and molecular ions with a 2 Tli/ 2 ground state (rep- 
resented here by group 15 fluorides). The values of Kqp 
for all the systems under study are presented in Table |H} 
together with the DHF and the DFT Wa constants. As 
the final recommended value for the Wa parameter we 
take an average of W a (DHF)Kcp and Wa{UFT)K C p- 
The estimate of the accuracy in our previous work [23] 
has shown that it is about 15% for molecules in the 2 E 1 / 2 
electronic state and 20-30% for the 2 Hi/ 2 state. 

The magnitude of Wa in the 2 £1/2 electronic state is 
expected to scale as Z 2 R\y [H], where Rw is the relativis- 
tic parameter, 



Rw 



2 7 + l 



(IB 



7 = [1 - (Za) 2 ] 



2Zr A 1 / 3 
211/2 



2-27 



[T(2 7 + 1)] 



(10) 



-15 



111, 



In Eq. (10 1, ap is the Bohr radius, r = 1.2 x 10 
and a is the fine-structure constant. The Rw parameters 
are shown in Table [TT| for each of the metal atoms. In 
Fig. 1 we plot log J j^^ • ) as a function of log(Z) for 
both groups of dimers. For group- 13 fluorides the scaling 
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FIG. 1. (color online) Log-log plot illustrating the scaling 
of Wa with relativistic parameter Rw and atomic number Z 
for group-13 and -15 singly-ionized fluorides. 

is, indeed, Z 2 . In the case of group 15 fluorides, however, 
the ground state is 2 Ui/2, for which the Wa parameter 



vanishes in the non-relativistic limit, since in this limit it 
does not contain the s-wave electronic orbital and can not 
provide the matrix element (si/2| Q: p( r )|Pi/2)- The effect 
appears due to the mixing of 2 S 1 / 2 and 2 n 1 / 2 electronic 
states by the spin-orbit interaction, and gives an extra 
factor of Z 2 a 2 in the zJ-dependence of Wa, as seen in 
FigB 

To summarize, here we have performed calculations 
of the P-odd interaction constant Wa in singly-ionized 
group-13 and group-15 fluorides as well as other select 
singly-ionized diatomic systems. To the best of our 
knowledge, this is the first investigation of nuclear spin- 
dependent parity violation effects in molecular ions and 
it is motivated by progress in the cooling and trapping 
of such systems. 
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